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■ Abstract: The plethora of scalar fields participating in the formulation of a softly bro- 

. ken supersymmetric theory can threat the stability of the standard vacuum. The generic 

, situation is twofold. Directions in scalar field space may exist along which the potential 

\ becomes unbounded from below or local minima deeper than the standard one are likely to 

\ appear where charge and/or color symmetries are broken. An investigation of these mat- 

I I ters requires a thorough study of the tree level effective potential as well as its radiative 

corrections. However, the presence of many different mass scales in such a generic super- 
' symmetric theory needs an appropriate renormalization improved treatment. Combining 

the decoupling theorem with a conveniently chosen renormalization scale for each field con- 
^ . figuration, the well known multimass scale problem is circumvented. In this context, the 

^ \ ordinary universal soft-parameters case of the minimal supersymmetric model as well as a 

non-universal case of a brane world minimal supergravity model are examined closely for 
unphysical vacuua. 
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1. Introduction 

Supersymmetric (SUSY) theories are the best motivated extensions of the Standard Model 
(SM) of the electroweak and strong interactions. They provide an elegant way to stabihze 
the huge hierarchy between the Grand Unification (GUT) and the Fermi scale, and predict 
a variety of new matter states (sparticles) creating a natural framework to cancel the 
quadratic divergences of the radiative corrections to the Higgs boson mass. One can classify 
these models by the mechanism for communicating SUSY breaking from a hidden sector 
to the observable sector. Possibilities include gravity mediated SUSY breaking (SUGRA) 
Q , gauge mediated SUSY breaking, Q and anomaly mediated SUSY breaking . 

The so-called minimal supergravity (mSUGRA) model has traditionally been the most 
popular choice for phenomenological SUSY analyses. In mSUGRA, it is assumed that its 
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most economical low energy realization, the so called Minimal Supersymmetric Standard 
Model (MSSM) |], is valid from the weak scale all the way up to the GUT scale Mx ^ 
2 X 10^^ GeV, where this choice is usually taken due to apparent gauge coupling unification. 
In the MSSM one assumes a minimal gauge group, i.e. SU(3)c x SU(2)l x U(1)y, and a 
minimal particle content, i.e. three generations of fermions (no right handed neutrinos 
included) and their spin zero partners as well as two Higgs doublet superfields to break 
the electroweak symmetry. In order to explicitly break supersymmetry while preventing 
the reappearance of quadratic divergences, a collection of soft terms Q is added to the 
Lagrangian: mass terms for the gauginos (^1/2)1 mass terms for the scalar fermions (mo), 
mass and bilinear terms for the Higgs bosons {Bq), and trilinear couplings (^0) between 
sfermions and Higgs bosons. In the general case, that is if one allows for intergenerational 
mixing and complex phases, the soft SUSY breaking terms will introduce a huge number 
of unknown parameters. This feature makes any phenomenological analysis in the general 
MSSM a daunting task. In this context, severe phenomenological constraints are imposed 
on the parameter space by flavor changing neutral currents (FCNC) and CP violation Q 
as well as unphysical vacuua. 

To reduce the number of free parameters, one needs an explanation of how super- 
symmetry is broken. In mSUGRA, supersymmetry is spontaneously broken via a hidden 
sector field vacuum expectation value (VEV), and the SUSY breaking is communicated to 
the visible sector via gravitational interactions. For a flat Kahler metric and common 
gauge kinetic functions /ab, this leads to universal values for mo, rni/2 ^Oi Bq at 
the GUT scale Mx- This assumption of universality in the scalar sector leads to the phe- 
nomenologically required suppression of flavor violating processes that are supersymmetric 
in origin. However, there is no known physical principle which gives rise to the desired 
form of Gj and Jab] indeed, for general forms of and Jab, non- universal masses are 
expected. Hence, the universality assumption is regarded as being entirely motivated by 
the phenomenological need for suppression of flavor violating processes in the MSSM. 

In recent years the branes, which are typical in models with extra dimensions have 
been found to fit naturally with the idea of breaking supersymmetry via a hidden sector. 
The basic five-dimensional setup of brane world models is that of four-dimensional hy- 
persurfaces (branes) hosting familiar gauge and matter fields, which are embedded in a 
higher dimensional ambient space, the bulk, populated by gravitational and gauge-neutral 
fields. The bulk degrees of freedom couple to the fields living on branes through various 
types of interactions, but the effects of SUSY-breaking will be communicated gravitation- 
ally to the observable world through the five-dimensional interior bulk. This scenario 
thereby places the question of SUSY-breaking in an entirely new geometric context. 

Brane world SUSY breaking generally gives rise to tree-level soft scalar masses for 
visible sector squark and slepton fields. These masses are generally not universal. Without 
additional assumptions about fiavor, these non-universal scalar mass matrices are not nec- 
essarily aligned with the quark and lepton mass matrices, and dangerous sfiavor violation 
can occur. A model to overcome this flavor problem has been recently proposed. The 
authors of Ref. |^ have suggested a gravity mediated SUSY breaking in 5D spacetime with 
two 4D branes Bl and B2 separated in the extra dimension. The SUSY breaking effects 
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from the hidden sectors locahzed on the branes are transmitted through gravitational cou- 
phngs both to the visible sector fields located on the same brane and to bulk gauginos. 
Thus, if the SUSY breaking scales are mi and ?Ti2 at two separated branes Bl and B2 
respectively, the localized chiral multiplets at Bl (B2) could get soft SUSY breaking effects 
of order mi (7712). Assuming that the SUSY breaking scales at Bl and B2 are mi (~ 10-20 
TeV) > m2 (~ 1 TeV), the superpartners of the first two generations can be easily made 
sufficiently heavy to overcome the SUSY flavor problem by locating them on Bl. In order 
to keep the radiative corrections to the Higgs masses under control, the third family must 
reside on B2. The two Higgs multiplets should reside in the bulk, so that the first two 
and the third generations of the quarks and leptons can couple to them at Bl and B2, 
respectively. 

Beyond the SUSY flavor requirements, an additional constraint on mSUGRA param- 
eters can be obtained by demanding that the global minimum of the scalar potential is 
indeed the minimum that leads to appropriate electroweak symmetry breaking |1C]. From 
a theoretical point of view, the plethora of new scalar fields which are introduced in MSSM 
may lead to many possible directions in field space where field configurations could be 
developed deeper than the standard minimum. The generic situation [ pH] , |l^, 14, 15] is 
that the scalar potential can receive negative contributions from quadratic or cubic term in 
the supersymmetry breaking Lagrangian if the usually dominated quartic D and F terms 
are suppressed. A systematic classification of all possible dangerous directions in the scalar 
field space that can potentially lead to undesirable minima has been done by the authors 
of Ref. |T^. These directions have been categorized either as field directions that are un- 
bounded from below (UFB) or as directions that lead to charge or color breaking minima 
(CCB). 

The purpose of this work is to extend and complement the work done on this sub- 
ject by analyzing the impacts the one loop radiative corrections have on these unphysical 
configurations for the models described above. Yet, the effective potential in which the 
vacuum structure is encoded, is a poorly known object beyond the tree-level approxima- 
tion. One reason for this is the dependence of its loop corrections upon the very many 
different mass scales present in these models, so that a renormalization group analysis 
becomes rather tricky. In general, when one deals with a system possessing a large mass 
scale Qm, compared with the scale at which one discusses physics, large logarithms 
such as ln((5A//<5/i) always appear which affect the convergence of the loop expansion. In 
this situation, one considers resuming the perturbation series by using the renormalization 
group equation (RGE). Nonetheless, in many realistic applications one often has to deal 
with an additional mass scale Qm with the hierarchy Qf^ <C Qm Qm- For instance, 
one can regard Qfj,,Qm,QM as the weak, super symmetry-breaking and unification scales 



respectively. When we study such a system, we face the problem of multimass scales [16|: 
There appear several types of logarithms ln{QM/Qf_i) and ln{Qm/Qf_i), while we are able 
to sum up just a single logarithm by means of the RGE. 

In the present work, trying to circumvent this problem, a generalized improving pro- 
cedure based on Refs [17, 18 1 is applied to the relevant effective potential. The main idea 
of the method is to make use of the decoupling theorem [19| combined with a conveniently 
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chosen renormahzation scale Q* for each field configuration. By this procedure, it is made 
possible to treat essentially a single log factor at any renormahzation scale, since all the 
heavy particles (heavier than that scale) decouple and all the light particles (lighter than 
that scale) yield more or less the same log factors. 

The rest of this paper is organized as follows. After setting our notation and conven- 
tions in Sec. |2|, an estimation of the physical vacuum is attempted in various effective low 
energy SUSY models in Sec. ^ Next, section ^ describes briefly the most dangerous direc- 
tions that could lead to unphysical situations and Sec. ^ is an outline of our calculational 
scheme. In Sec. ^ we give a brief summary of our results and section summarizes our 
conclusions. Finally, detailed formulae for the various mass matrix elements involved in 
one loop expressions of the effective potential are presented in Appendices. 

2. The Physical Set-up 

Let us briefly review some of the basic ingredients required for our analysis. A globally 
supersymmetric and SU{'i)c ® SU{2)l ® U{1)y gauge invariant Lagrangian with minimal 
content can be constructed from the usual i?-symmetry conserving superpotential^ 

W = Y,VE'n\eij + YdQ^''b'in\eij + Y^Q^^U^nieij + ^,n{nie^J (2.1) 

where = (n d), U'^, D'^ are the Quark Superfields, V = (z> £), E'^ are the Lepton 
Superfields and ?Yf = {Hi hi), = {h2 H2) are the Higgs Superfields. Generally the 
Yukawa matrices Y^, Yd, Y^ and the parameter ^ are complex. 

In order to explicitly break supersymmetry as required by experiment, while keeping 
quadratic divergences suppressed, a collection of soft terms is added to the Lagrangian. 
These include: mass terms for all scalar fermions and gauginos, bilinear terms for the 
Higgs bosons and trilinear interactions between sfermions and Higgses. Consequently, the 
most general soft SUSY breaking Lagrangian with real mass terms in this minimal scheme 

is2 

C-soFT = -\mi{BB) - ]^M2{W^W^) - ^M3(G«G^) 

+huQ^''u;^nyij + H.c.) - {Bfininieij + h.c.) (2.2) 

Here 7ii, 7i2 are the ordinary Higgs boson doublets, Q, D^, U^, L, E'^ are the squark- 
slepton scalar fields and h = Y A, where A is a 3 x 3 matrix containing the "soft" trilinear 
scalar couplings. All extra soft parameters except masses are generally complex. 

Altogether one would then need more than 100 real parameters to describe SUSY 
breaking in full generality. Clearly, some simplifying assumptions are necessary if we want 
to make comprehensive scans in parameter space. Specifically, we shall work with real 

are SU{2) indices, a is a color index (family indices are suppressed). Also ei2 ~ +1. 
^B, W, G stand for gauginos (Weyl spinors) and F, R axe SU(2), SU{3) group indices respectively. 
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and diagonal Yukawa couplings where all non-zero entries are positive. Also B and 
all trilinear soft couplings A are assumed to be real. In our analysis we shall also keep 
Yukawas and trilinear soft couplings from light families, since their contributions to the 
one-loop effective potential (our main objective) are not always negligible for an arbitrary 
field configuration. For clarity reasons, we simplify the notation using yj = Y^, = 
Yt, YJ^ =Y(., Y^ = Yh and similarly for the trilinear couplings. 

A dramatic decrease of the unknown parameters arises when the minimal model just 
presented is embedded into various supergravity scenarios. The simplest case of the MSSM 
results from coupling with = 1 Supergravity. This leads to the following "universal" 
scenario at a very large scale Mx — 2 x 10^^ GeV with: (1) Common gaugino mass : mi/2j 
(2) Common scalar mass : mo, (3) Common trilinear scalar coupling : Aq. This reduces the 
number of free parameters describing SUSY breaking to just four : The gaugino mass 1711/2^ 
the scalar mass mg, the trilinear and bilinear soft breaking parameters ^0 ^-nd B, which 
is conventionally trade with the Higgs VEV ratio tan/3 = {H2)/{Hf). We also assume 
unification of the gauge couplings at the scale Mx, while no specific relation is assumed 
for the Yukawa couplings there. 

However, more complicated alternatives also exist which give definite predictions for 
scalar, gaugino masses and soft trilinear terms appearing in Eq. (p.2|). As we have men- 
tioned in the introduction, one such possibility is the recently proposed 5D supergravity 
model of Ref. . In this model the chiral matter content resides on two 4D branes sepa- 
rated in the extra dimension. Higgs bosons as well as gravity and gauge multiplets reside in 
the bulk. It is then shown that the SUSY breaking effects from the hidden sectors localized 
on the branes are transmitted to the visible sector fields located on the same brane and 
bulk gauginos through gravitational couplings. Trying to overcome the SUSY flavor prob- 
lem, two SUSY breaking scales mi,m2 are introduced at the separated branes Bi and B2, 
respectively. Then the localized chiral multiplets at i?i(i?2) could get soft SUSY breaking 
effects of order mi (7712). Assuming that the first two generations reside at Bi and the third 
one is at B2, the following non universal soft terms are derived 
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-{Sa + 5,2) + -r^5i3 (2.3) 



M = ^{S^, + 6^2) + ^S^s (2.4) 

where i runs on three generations. Mass terms for the gauginos can be generated at tree 
level in a D = 5 off-shell SUGRA formalism |20|. If SUSY is spontaneously broken at 



brane Bi, heavy gaugino masses given by mi/2 — mi {M^/Mp), where M^, is a 5D SUGRA 
fundamental scale and Mp is the reduced 4D Planck scale, could give rise to large radiative 
corrections to the Higgs mass and spoil the naturalness solution. However, a not so "large" 
extra dimension, for example of the order of L^^ ~ 10^^ GeV (or Af* ~ lOL^^ ~ 10^^ GeV) 
suppresses the gaugino masses to 1 TeV scale. Note that the cutoff scale Ac introduced to 
obtain an effective 4D theory from the 5D one, is of the order of the extra dimension. Hence, 
the effective theory below Ac ~ Mx resembles the familiar gravity mediated scenario of 
MSSM. 



- 5 - 



In order to discuss the physical implications of these models at low-energies, we have 
to determine the evolution of the various "couplings" at experimentally accessible energies. 
For that reason, we use two loop RGBs with all thresholds of the various particles included 
|21|, E2|. 



3. Estimating the physical vacuum at low energies 

Regardless of the physical scenarios describing physics at high energies, low energy world 
could be best described by a softly broken super symmetric theory. In such a framework, 
the elegant ideas of supersymmetry, unification and radiative symmetry breaking coexist. 
The minimal super symmetric extension of standard model (MSSM) incorporates all the 
above. Because of its minimal content and the radiatively induced symmetry breaking, it 
is the most predictive of analogous theories. Let us denote by cf) the scalar field collection of 
MSSM. Since MSSM is expected to be a stable theory, its full effective potential Vfuii (scale 
independent) must have a global minimum (physical vacuum). At the physical vacuum, the 
only scalar fields that should acquire non-trivial expectation values (VEV) would be the 
neutral Higgs components. Simply speaking, the VEV configuration there should read:^ 
{{Hi), (Hi); 0). 

We know that MSSM depends on a minimal set of free parameters {Aq, mo, ^-1/2; 
tan/3). On general grounds one expects that, for some values of them, \fuii might be 
unbounded from below (UFB) (no physical vacuum). For others, a global minimum may 
exist where additional scalar fields acquire non-trivial VEVs breaking in this way at least 
one of the other gauge symmetries {U{1), SU{3)), a fact that clear contradicts experiment. 
How should we protect the theory from cases like that? The most straightforward way is 
direct minimization of Vfuu for each "point" (^O; 1^0 ^ t^-n P) and exclusion of the cases 
that give non-physical results. 

Let us employ this procedure for a physically stable "point" where only SU{2) is 
broken. To compute the VEVs we group the scalar fields into two subsets: = {Hf;xk} 
where i = 1,2 and k runs on the rest scalars. Applying the stationary conditions we get: 



and 



= (3.1) 



From a more practical point of view, one must construct an approximation of these abstract 



conjectures using loop expansion Indeed, at one-loop level the effective potential p4| 



becomes scale dependent and we can write 

Vfuii ~ Vi^ioop{Hl Hi, x; Q) = n' + y w + (3.2) 



where Q is the renormalization scale and is a field independent quantity |^5[ that should 
be added to preserve scale invariance at one-loop level. On the other hand, V^^^ and V^^^ 
are the classical and one-loop corrections to the scalar potential respectively. 



^Since \juu is 5(7(2) symmetric, we can choose either (-ffi) or {H2) to be a real non-trivial number. 
Then the reality of the other should arise from the minimization conditions. 
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Since V^'^'^ is at least a second degree polynomial in Xfcj its first derivatives with respect 
to Xk vanish at physical vacuum and the 1-loop approximation to Eq. will become 



dVi^loop{HlHlO;Q) 



dH? 



and 



dXk 



(3.3) 



In order to proceed further, we consider the "reduced" function Vi^ioop{Hi, H2,0;Q) and 
expand around the exact one-loop solution {vi,V2), {Hf) = vi + 6vi, {H2) = V2 + 6v2 
where 5vi are of higher-loop order. We get 



dVi.ioop{H^,HlO;Q) 



dH? 







(3.4) 



Solving this system of equations, we get the one-loop approximation of the so called physical 
vacuum. As an example, we will present the well known expressions for the physical vacuum 
at tree level. The associated "reduced" function is 

V^'\HlHlO-Q) = ml\H^\^ + ml\Hl\^ + 2mliH',Hi) + ^±^{\H^\^ - l^/^T)' (3-5) 



and its value at the minimum is ||13| 



V 



(0) 



ph 



[ml — m\ + {m\ + 7712) cos 2/3] ' 



(3.6) 



where tan (3 stands for the VEV's ratio at tree level. 

We stress here that the "reduced" function Vi^ioop{H^i 0; Q) appearing in ( |3.4D can 
be used solely as far as the computation of vi,V2 is concerned. It is completely unsuitable 
for drawing any conclusions about the kind of l^uii's stationary points. Nevertheless, these 
"reduced" functions can help us to study more general cases. Practically, one selects a 
"direction" in the field space"^ whose "reduced" function collects the deepest values of Vfuu, 
at the approximation level we work, and excludes its trace from the model's parameter 
space whenever it leads to a non-physical picture. 



4. Classification of dangerous directions in a nutshell 

Scalar potential in SUSY models receives contributions from three sources: D-terms, F- 
terms and soft-breaking terms. L'-terms provide a quartic contribution V ~ Xcj)^ with 
A > 0. The special case A = can only occur along special directions in field space known 
as "D-flat directions". Along such a flat direction the quartic potential takes the form 
V ~ (0f — (/>2)^ ^ as |(/>i| |(/)2|, which may render running off to infinity. Whether 
or not this occurs, will depend on the magnitude of the other contributions (F-terms, soft 
terms). 

F-terms, in turn, due to their supersymmetric nature contribute quadratic, cubic, and 
quartic terms to the potential. This does not mean, however, that the minimum of the 

*0f course, now all the scalar fields are allowed to participate 
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potential lies at the origin. Rather, it means that directions in field space with non-zero 
quartic F contributions will be well-behaved far away from the origin. Still, there will be a 
subset of the D-flat directions which are also F-Hat and whose behavior will be completely 
controlled by soft-breaking terms. 

The soft mass contributions, on the other hand, are problematic, as they can have 
either sign. Because they contribute to only the quadratic and cubic pieces of V, one can 
analyze their structure most readily along flat directions in which the quartic pieces all 
vanish. Then, one finds two distinct types of problems which may arise: potentials that 
are unbounded from below (UFB), or potentials which break charge and/or color (CCB) at 
their minima. Following Ref . ^] , we will present a tree level based classification of the 
various dangerous directions in various low energy supersymmetric models with universal 
or non-universal soft terms. 

4.1 UFB cookbook 

In order to present the most dangerous of the UFB directions, it is useful to notice some 
general properties:^ Soft trilinear scalar terms AYFHF^ cannot infiuence dramatically a 
UFB direction, since for large enough values of the fields the relevant F-terms dominate. 
The only negative contributions in V^'^'^ that can affect a UFB direction are:^ m2|//2p, 
— 2|m3||ifii72|- Since these terms are quadratic, along a UFB direction all the quartic 
positive terms coming from F- and D-terms must be vanishing or kept under control. 

• UFB-1 The first possibility is to play with just Hi, H2. This is the ordinary example of 
a UFB direction found in the neutral Higgs potential along the direction \Hi\ = \H2\. As 
with all UFB potentials, there is no quartic contribution to the potential along this direction 
(a cubic piece is also absent here). In the present framework, since we are mainly interested 
in the effect of one-loop corrections, we will also allow the charged Higgs components to 
participate. This is calculationally tractable only for this case contrary to the subsequent 
ones. Our choice for the "reduced" function will be'^ 

VUFB^I =V^'^\Hi,hi,H2,h2) + V^^\Hi,hi,H2,h2), (4.1) 
= ml {\Hi\'^ + |/ii|2) + ml (|F2|2 + |/i2p) + 2ml {H1H2 - /ii/i2) 
+ I \H1h2 + h,H2\' + + \h,f - \H2\' - \h2\')' (4.2) 

Z o 

and necessary ingredients are presented in Appendix |C|. 

• UFB-3 One other possibility is to take Hi = 0, H2 7^ 0. We shall try to control 

I 1 2 ~ 

the term \nH2\ by introducing some squarks or sleptons. Notice that a nontrivial U^^ 
(specific r, a) will produce dominant i^-terms, so it must be excluded. Using the fact that 
= { iL d) and L?" = { v £ ), the remaining i^-terms are 



y,r<F2r + 



2 



+ 



fiH2 + YJlrE^ + Y^dP^D. 



2 



^In what follows H1^2 ~ Hi^2, — hi, H2 = h2. In addition r,s,t stands for family indices and i,j\a 
are SU {2); SU {?>) indices respectively. 



^Throughout, we use the convention m% = < 0. 
With no harm of generality, we will assume real fields. 
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If we impose v!^ = 0, then the non-trivial terms will become 



fiH2 + Y^irE^ + Yjd'^D',. 



Further suppression can be achieved in two ways: 

a) E'^ = 0, and Df.^ ^ and 7^ for a specific value of r, a. 

In this case we do not eliminate quartic terms, but we try to neutralize their effect 
using additional scalars. It is not difficult to see that the strongest constraint arise 
for r = 3 (a whatever). An intuitive argument in favor of this is considering the 
soft squark contribution m'^^^\D^\'^ to V^^\ As we approach lower values of V^^\ the 
relevant quartic F-term deactivation leads to \D'^\'^ —>■ \fiH2\/Yii. So, sbottom gives 
the less positive contribution. Control of D-terms also necessitates the presence of 
some slepton doublet provided that it does not introduce additional F-terms. Hi or 
charged Higgs hi,h2 cannot play this role due to non-trivial quartic terms. Thus, 
the simplest choice which gives a minimal positive soft contribution and suppresses 
the quartic D-term \H2£r\^ is LJ = { Vr ) 7^ for some specific r. Our "reduced" 
function will then be 



VuFB-^a = V^'>\H2,d3,Dlh)+V^'\H2,d3,Dl93), (a = 1) 



(4.3) 



= {ml-^^')\H2\' + ml\Ds\' + ml^\ds\' + m%.\DI\^ 



+ 



fiH2 + Ybd^D'^ 



j2 1 2 



3I i\ 



-\D 



c|2 



\D 



c|2 



(4.4) 



b) D'^^ = 0, and / and D^- = and F^ 7^ for some specific r. 

Similarly as in the case (a), the strongest constraint comes for r = 3. Further control 
of D-terms requires a slepton doublet from another family {h2 is inappropriate since 
it introduces quartic terms). So we take = { Us ) 7^ and F^ = for some 
specific s r. Choosing s = 2, we get the following "reduced" function 



V, 



UFB-3b 



yW(i/2, 4, ^2) + V^'HH2,£3, El D2), 



(4.5) 



= {ml-^^'')\H2\''+mlfU^+ml}u2?+ml^^\^ 



+ 



flH2 + YrlsE'. 

2 



+ ^ 



' ,~, |2 , Iff |2 



+ — I |Z^2 



1^3 



^2r + |^3| 
2 



2\E'. 



•c\2 



(4.6) 



while pieces are shown in Appendix O 
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4.2 CCB cookbook 



CCB most readily occurs along directions which are £)-flat, though not necessarily F- 
flat (the (j)'^ contributions to the potential are suppressed by Yukawas). In order to gain 
intuition about CCB, let us present some general comments: The most dangerous CCB 
directions involve only one particular trilinear coupling of one generation. Two or more 
trilinear couplings with different values of Y can not interfere constructively in the same 
region of field space to deepen the potential. The CCB directions we explore are not F-flat, 
so the most restrictive are those with the smallest non-vanishing F-terms. Since they are 
proportional to the respective Yukawa coupling, we expect cases with the lightest Y to be 
more restrictive. The most dangerous CCB directions are: 

• CCB-Ea Our goal here is to keep only the "light" trilinear coupling. For that 
reason we assume Hi ^ 0, ^ 0, LJ = { i>j. £r ) ^ (for some specific r). The relevant 
F-terms are 



Y^DrE^f + \Y^£rE^\^ + l/^^il^ + \Y;HiE;^\'' + [Y^lrHil 



'?c\2 



Obviously further suppression is possible so we choose i^r = (the case ir = kills the 
Ae trilinear term) and introduce some squarks (sleptons will add a new A) to control the 
pure Higgs F-term. We can take Qf^ = { Ug dg ) ^ 0, and ^ (for some specific s, a). 
Then the relevant F-term becomes 



and the minimal choice is for dg = 0. Using a similar intuitive argument as in UFB case, we 
conclude that Yt gives the smaller positive F-term. Choosing r = 1 (lighter Yukawa), s = 3, 
a = 1 makes further control of D-terms unnecessary. Putting altogether, the "reduced" 
function becomes 



VccB-Ea = V^'^HHi,£i,E'i,U3,Ul) + V^'\Hi,£i,E'l,us,U^), 



(4.7) 



F(o) = (m?-/x2)|i/i|2 + m|,|Fj|2 + m|j^~i|2 + m|jti3|' + m|,|C/||2 



+ 2YeAj:iE'iHi 

■,2 



+ ^ ( 1^1, 

^2 



2\E' 



+ 



1(141 



\U3 



YjiEll^ + lYeElHif + \YeiiHi\ 

2 



(4.8) 



• CCB-Eb The other way to keep A^, is by two Higgses Hi 7^ 0, i?2 7^ 0. Now one should 
allow the same sleptons E^ ^ 0, LJ = { i>r ir ) ^ (for some specific r). However, here we 
cannot use squarks or sleptons to compensate the Hi F-term, because this will introduce 
an additional trilinear coupling. The relevant F-terms will be 



|y;p.F,f + 



+ i^iHif + \y;e^.Hi\^ + lYjirHif. 
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By taking = 0, we ensure minimal i^'-terms without kilhng the triUnear contribution. 
In the Ughter case (r = 1), further control of D-terms is not necessary, so our "reduced" 
function will be 

VccB-Eb = V^'^\Hi,H2,li,E^) + V^^\Hi,H2ji,El), (4.9) 

= ml\Hi\^ + mil/fal^ + ^J^^iP + m\j,El\^ + 2mlHiH2 

+ 2YjiEl{fiH2 + AeHi) + \YjiEl\^ + \YeElHi\^ + \YeiiHi\^ 

+ ^ ( i^ii' + - - \H2fY + 4 ( \H2\' + - m'YiA.w) 



Expressions of V^^^ are presented in Appendix ^. 
5. Calculational Scheme 

Let us now describe the method used to probe the dangerous directions just presented. As 
is well known, one-loop effective potential depends on the eigenvalues of the tree level mass 
matrices. To be more accurate, these corrections are derived under the assumption that 
the potential is resting at a tree level minimum (eigenvalues of dijV^'^^ positive). However, 
in a UFB case we need to define the one-loop corrections away from a classical minimum. 
Hence several eigenvalues of the scalar mass matrix may be negative leading to a complex 
valued function. These contributions are a signal that the sum of one-particle irreducible 
diagrams Vipi does not give the effective potential. Formally, in a case like that we must 
use the convex envelope of Vipi, which in the vicinity of a classical minimum matches the 
usual loop expansion. Trying to estimate the one loop corrections despite non-convexity of 
the effective potential, we adopt a moderate approach, namely V^^^ everywhere is given by 
the real part of the ordinary one-loop expressions [26, p^. Thus, in a mass- independent 



renormalization scheme (DR) V^^^ is given by (Q = Qe^^^) 



vW=k Y: y}'^ with y/^) = ^^(25, + l)C,AA,M^4(<^)ln^S^ (5.1) 

(Af2^0) 

where k = {IGi:'^)^^. Alp stands for the tree level mass eigenvalue of the p*^ particle, J\fp 
is the number of its helicity states and the associated spin, color degrees of freedom are 
denoted by Sp, Cp respectively. It is evident from ( |5.1D that in the case of a single mass 
scale a judicious choice of the renormalization scale (Q^ = M'^^cp)) eliminates all large logs. 
However, in the case of many different mass scales any renormalization scale will leave 
large logs remnants behind, so that we need higher loop corrections in the loop expansion 
to trust the results. 

The heart of the problem lies in the renormalization scheme we use (DR). For a mass 
independent scheme the decoupling of various mass states is not automatic and has to 
be incorporated. Hiding all heavy particle loop contributions in a redefinition of some 
lower energy parameters, we secure that all masses smaller than a given scale (decoupling 
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scale) behave as massless, while larger masses decouple and never generate problems in 
perturbation theory. Below a decoupling scale the theory is an effective field theory with 
new RGEs, while threshold effects take care of the matching between both theories at the 
boundary scale. 

A simple way to realize this scenario is to treat the thresholds as steps in the particle 
content of the RGE /3-functions |^2[ . Usually one integrates the RGEs from a superlarge 
scale Mx to any desirable value of Q- As we come down from Mx, as long as we are at 
scales larger than the heaviest particle threshold, we include contributions from all particles 
in the model. When we cross the heaviest particle threshold, we switch to a new effective 
field theory with the heaviest particle integrated out and of course a new (3. For field 
configurations in the low energy regime (< 300 GeV), the condition to determine the exact 
point of decoupling is simply = \w?{Q)\, where m?{Q) is the running soft parameter 
corresponding to the particle.^ Obviously, the step functions in RGEs will have the form 
Om = 0{Q'^ — \m?'{Q)\). Alternatively, for all other field configurations, the decoupling points 
are fixed by = \M'^{(f); Q)\, where is the field dependent mass eigenvalue of a particle. 
Analogously, the step functions in RGEs will now become Om = 0{Q'^- \M^{(j); Q)\), where 
(p are in principle all the relevant fields along the dangerous direction under consideration. 
Additionally, knowing that the minimum of hes at a non trivial field configuration, a 
suitable generalization of vacuum subtraction in Eq. ( |3.2D is given (at one-loop order) by 



m 



-n' = y(0)((i/i), {H2);0) + F(i)((i?i), {H2),0; Q) (5.2) 
Finally, we also replace the potential (^]^) by [^] 

V^^^ = kY^V^^'^6i where e^ = e{Q^ - \Mf{^-Q)\^ (5.3) 



One other issue is that of choosing the renormalization point. Following |18], we use 
a prescription that preserves the hierarchy (|T^*-''^| > l^^^^l) a-nd has a continuous field 
dependence 

Q* = 10'^(^)v/l^, $ • * = ^ 02 (5.4) 

i 

where x = log(\/$ • $/A) (A = 1 GeV makes the log argument dimensionless) , (pi are the 
scalar fields involved in the dangerous direction under consideration and our ansatz for 
iv{x) is given in Appendix ^. 

We stress here that not only the effective potential in Eq. (5.3), but also the (3 functions 
in RGEs do depend on the same step functions 



d\nQ 



/3^(0,(Q2_|m2(0;Q)|)) (5.5) 



Formally, we can integrate these implicitly field dependent RGEs to the appropriate renor- 
malization scale Q* and then construct the required effective potential. However, from 



a practical point of view using all the relevant scalar fields involved in Eq. (5^) requires 



*We use the factor e^'''' for compatibility with Eq. ( ^.l| ) 
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excessive computational time. The computation of the field dependent mass eigenvalues 
Mf{<j);Q) is very complicated due to higher dimensional mass matrices and has to be re- 
peated for each integration step taken internally by our numerical integrator. A possible 
way out is to restrict ourselves to the neutral Higgs subset of (/>'s solely in Eq. (5.5). In that 



case the mass eigenvalues are trivially given by solving algebraic equations, but the most 
important is that the computational time is reduced by several times. We stress here that 
this kind of approximation is applied only to the /3 functions in RGEs (Eq. ( |5.5D ) and not 
to the effective potential (Eq.(^^)) otherwise loop expansion is in danger. We have also 
checked the results of the two alternatives just described for some representative points of 
the parameter space and found no significant difference. So, finally in our investigation we 
have adopted the second way which is the quickest. 

To complete the picture, one also needs some "boundary scale" Qhigh where the starting 
values of the running parameters (couplings, fields) should be provided for the evolution at 
Q* . Notice that now, due to field dependent thresholds in /3 functions, the RGE evolution 
for Q < Qhigh depends on the field point we are. A convenient choice for Qhigh, besides 
the unification scale, is an intermediate scale higher than the largest field dependent mass 
eigenvalue at the current field point. Valid choices for Qhigh are Qhigh ^ 2.12000, where (poo 
stands for the upper bound order of magnitude of the allowed values for the scalar fields 



1 18]. Besides, we also need to know the values of the running parameters and fields there. 
For the fields the most plausible option is to take 4>{Qhigh) = fpo where (po, is the field point 
we examine. Since Qhigh is above all thresholds, the required values for the couplings there 
should not depend on the background fields and a reasonable choice is 

KiQhigh; M = KiQhigh; {Hi), {H2),o) (5.6) 

where the RHS is obtained by running the couplings from their known values at 
Mz when the potential is resting at its physical vacuum. Evolving this set of values 
{4>o, XaiQhigh'i {Hi), {H2),0)} from Qhigh to Q* using field dependent thresholds, the ef- 
fective potential at the current field point can be constructed. 

Minimization of the effective potential is performed numerically. Since analytic expres- 
sions for the derivatives are not available, we resort to methods that require only function 
evaluations. One such very efficient method is the downhill Simplex method |2^. A sim- 
plex (or polytope) in n dimensional Euclidean space is a construct with n + 1 vertices 
defining a volume element. For instance, in two dimensions the simplex is a triangle, in 
three dimensions is a tetrahedron and so on. Using a population of n -|- 1 points (simplex 
vertices), the algorithm brings the simplex in the area of a minimum and adapts it to 
the local geometry. The initial Simplex may be constructed from the current point (first 
vertex) by taking a single step along each of the n dimensions. 

6. Results 

Using the procedure outlined above and Merlin package |Q, we explored regions of MSSM 
parameter space for unphysical vacuua. We begin our discussion for the allowed parameter 
space in the {mo,mi/2) plane. We fix tan/3 to be 2 or 10, and take Aq = 0. Our search 
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Figure 1: Exclusion plots in the mo vs t7t,i/2 plane for the UFB-3b direction 

was performed in the range (0, 500) GeV for both axis and the grid was scanned with 25 
GeV resolution. In Fig. |l| the shaded regions to the left of the solid line display points 
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Figure 2: Exclusion plots in the mo vs Ao plane for the UFB-3b direction 
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where UFB-3b escapes were discovered. These parameter values with mo smaller than 
777-1/2 indicate that much of the mSUGRA parameter space associated with light sleptons 
is ruled out. It is also worth noting that higher values of tan/3 weaken these constraints. 

Although Fig. |l] is plotted for Aq = 0, a similar excluded region results for other 
choices of the parameter. An analogous plot in Fig. Q displays the effect of variation in 
the trilinear Aq parameter in combination with niQ for fixed mx/2 at 200 GeV. Again, as 
usual, we fixed tan/3 to be 2 or 10. Performing a numerical minimization along the UFB- 
3b direction, unphysical configurations were found only for initial values taken from the 
shaded regions. Similarly, we find out that the vacuum constraints disfavor light slepton 
cases. 

We have also examined a more general non-universal case coming from string models, 
which encompasses the special case where supersymmetry is broken in the dilaton sector. 
In the latter case one is led to GUT scale soft terms related by mi/2 = — = V^mQ. 
However, in our examination we treat tuq as a free parameter. We performed scans in the 
(mo,mi/2) plane with the usual values for tan/? (2 or 10). Numerical minimization showed 
that for small mo much of the parameter space is excluded by the UFB-3b constraint, as 
shown in Fig. ^. Furthermore, it is easy to see that the so-called dilaton dominated scenario 
corresponds to a straight line located entirely inside the forbidden region. 

A similar scan in the parameter space of the MSSM has also been performed along 
the dangerous direction CCB-Eb. We considered again the three different regions depicted 
in Figs ||-^. Using the simplex minimization procedure and the one-loop expressions in 
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Appendix we found no dangerous CCB minima. Tree level constraints described in 
Appendix ^ have also considered, as the numerical procedure scanned various field con- 
figurations and found that they are not violated for fields ~ A(,/Ye, ~ 10^/10~^ GeV 
~ 10® GeV. In great measure this is due to the large and positive contribution to the 
potential of the soft masses and especially m^. Let us stress here that our renormalization 
scale choice depends on the field point we are. If at some field point (pQ the renormal- 
ization scale Qq = Q*{4>q) does reverse the inequality of a tree level constraint, then the 
tree level potential constructed from "soft" parameters at the specific scale Qq will be 
deeper than the physical vacuum for fields ~ ^e(Qo)/^e(Qo)- In other words, tree 
level constraints signal an unphysical situation only if they are violated for fields of the 
order ~ A,{Q*) /Ye{Q*). 




Figure 4: Exclusion plots for the mi vs TO2 plane in a brane world model. 

We next focus our attention on the non-universal brane model and perform a scan for 
unphysical vacuua. In that case, due to non-universality, sparticles may generate too large 
flavor-changing (FC) effects but the specific brane model offers a mechanism for suppressing 
these dangerous processes. In this scenario, pushing the SUSY breaking scale at the first 
brane high enough (mi ~ 10 — 20 TeV), the masses of the first and second generations 
scalars will become sufficiently heavy to avoid the flavor problem. On the other hand, 
introducing a second SUSY breaking scale (m2) attached to the second brane constrains 
the heavy scalars to be not much heavier than 1 TeV or so to preserve the gauge hierarchy. 
We have performed scans in the {mi, 1712) plane for the values described above. Our results 
for tan (3 = 2, 10 are shown in Fig. ^ where one can also see the regions excluded by 
phenomenological constraints^ or/and by unphysical vacuua driven radiatively by negative 
scalar mass eigenvalues. Note that in the case tan/3 = 10, fi > the parameter space is 
strongly restricted allowing only a very small "window" for phenomenology. 

^By this we mean either no electroweak breaking or no mass for the gluino or violation of the experimental 
bounds for superpartners masses 
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Unfortunately, an analytical treatment of these issues is tedious and leads to very 
unwieldy expressions |31|. Relaxing our preciseness requirements, we can intuitively realize 
the strength of a dangerous direction by resorting to the familiar tree level expressions. 
Constraints based on tree level logic have been thoroughly studied in the past [12, 13, 15| 
and our intention here is not to reproduce these analyses, but merely to use them as an 
alternative confirmation since our renormalized effective potential respects perturbation 
series hierarchy. Indeed, as one can see in Appendix the global minimum of Eq. ( [4.2| ) 
(UFB-1 case) is not deeper than the physical vacuum, a fact that has been also numerically 
justified at one-loop level. 

One way to test reliability of the effective potential is by checking scale invariance 
1 27, In our case the chosen scale^'^ can not vary arbitrarily, since this endangers 

convergence of loop expansion. For very large values of scale 6 steps in V^^^ allow all 
heavy masses and loop corrections become huge. We have tried several scale prescriptions 
(various oj choices) and found that only the location of a UFB-3b escape is sensitive to 
scale variations. We believe that this sensitivity is due to cancellations that take place 
along a UFB-3b escape and higher loop effects. Since our expressions are only one-loop 
scale invariant, the higher loop difference is amplified by these cancellations and appears 
as a small "rotation" of the escape trajectory. The important point is that the physical 
picture remains unchanged: whatever lo{x) one chooses (within perturbativity constraints) 
cannot render an unstable potential stable and vice versa. 

In the present work we mainly focus on an explicit treatment of the effect one loop 
radiative corrections may have on various UFB and CCB dangerous directions in super- 
gravity models. Due to the intrinsic complexity of this subject, as well as the theoretical 
problem of multiple mass scales these issues have not been straightly touched upon in the 
past. Here for the first time, implementing a sophisticated threshold technique |17, 18 1 
outlined in Sec. |5|, we were able both to overcome the multiscale problem and deal ex- 
plicitly with the radiative corrections. We believe that this is a more complete approach 
leading to more reliable results and corresponding plots. Practically, our improved results 
are in broad agreement with previous investigations ||l^, |l5|. However, a more detailed 
comparison is highly non-trivial due to the variety of methods implemented by the various 
authors. 



7. Conclusions 

Minimal SUGRA scenarios provide a well motivated and phenomenologically viable frame- 
work of how weak scale supersymmetry might occur. Unfortunately, one has to pay the 
price of introducing extra degrees of freedom and new parameters leading to the appear- 
ance of new sources of flavor changing and CP violating processes, large number of new 
parameters and complicated scalar sector of the theory. On the other hand, the way low 
energy is related to the fundamental theory may shed some light to the solution of these 
problems. 



Note that our renormalization point is field dependent rendering Q,' in Eq. (5.2) a non trivial subtraction 



-17- 



In this paper, trying to constrain the parameter space of some mSUGRA models as 
described in introduction, we follow a strategy of finding regions in the general scalar field 
configuration space which have inferior values than the physical vacuum. The argument is 
that, if this configuration is not a minimum, an unbounded from below direction or a global 
minimum with charge and/or color broken surely exists somewhere and the associated point 
in parameter space should be excluded from consideration. 

We have analyzed the relevant potentials employing the full one-loop radiative cor- 
rections in the calculation. Because of the various mass scales present in these models, 
renormalization group improvement of the potential has ambiguities and should be care- 
fully treated. Implementing the decoupling theorem in a manner proposed by the authors 
of Ref. we treat the various particle thresholds as steps in the /^-functions as well 

as in the one loop corrections of the scalar potential. We stress here the role played by 
the renormalization scale choice, as given in Sec. It should be wisely chosen in order to 
eliminate heavy particles whose participation puts in danger the convergence of the loop 
expansion {i.e. < 

Employing the framework just stated and using numerical minimization procedure, we 
have performed an analysis of how the most dangerous directions put restrictions on the 
whole parameter space of the various models described. In both the MSSM case and the 
brane world scenario considered, these constraints turn out to be very strong, producing 
important bounds not only on the value of A (soft trilinear coupling), but also on the values 
of mi/2 (gaugino masses) and niQ (scalar masses). Our analysis is summarized in Figs 1-4. 
As a general trend, the smaller the value of ttt-q the more restrictive the constraints become. 

Finally, we note that the more or less strict constraints derived here from non-standard 
vacuua are avoided if we adopt the idea (S^ that we may indeed exist in a false vacuum 
and the tunneling rate from our present vacuum to a non-standard one might be small 
relative to the age of the universe. However, it is difficult to realize the circumstances 
under which the idea that we may live in a false, metastable vacuum could be reconciled 
with the existence of a small positive dark energy in the universe, either in the form of a 
constant vacuum energy/cosmological constant or in the form of a new scalar quintessence 
field, as the recent discoveries suggest [^ . 
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A. Analytic treatment of UFB-1 at tree level 

In this appendix we compute the global minimum of the tree level effective potential along 
the UFB-1 direction. Assuming for simplicity real fields, the solution is based on the new 
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"polar" variables {Ri,R2,^p,0) where 

Hi = Ricosip hi = Risimp H2 = R2Cos0 h2 = R2sm9 (A.l) 

Of course i?i > 0, -R2 > for an invertible transformation, thus the trivial case should be 
examined separately. If we take both or at least one R trivial, we also find a trivial value 
for the potential extremum which is clearly above the physical vacuum shown in Eq. 
Hence it is sufficient to minimize the function 



2 2 2 

mjRj+mlRl+2mlRiR2Cos{ip+e) + ^RlRlsm'^{^+e) + ^ '^"^ ^^^^ 



-{Rl-Rlf (A.2) 



Solution of the minimization conditions leads to the following cases: 



(I) R1R2 cos(v3 + 9) = 2mllgl, ^^^^^ = and dR^V^^^ = with i?2 ^ 
Solving with respect to Ri,R2 we get -Rf = {m\ — m^jg^ 



and 



i?2 = ("^1 — ^2)10^ — {m\ + mD/g'^. However, stability along Ri = ±i?2, 6 = cf) = Q 
direction dictates positiveness of the m\ + m| combination (m^ + > 2|m3|) so 
one of i?f , R2 will be negative. Thus, no stationary solution exists or the potential is 
unbounded from below. 

(II) sin((^ + ^) = 0, dn^V^^^ = and Br^V^^'^ = with i?i,i?2 / 

In our conventions we always have m| < 0, so smaller values for F^'^^ are obtained 



for cos((/9 + 0) > 0. So we should take here cos((^ + 6) 
the extremum conditions will become 



1 (i.e. 9 = —<p). The rest of 



mfi?i +m|i?2 + g^RiiRj 
mlR2 + mlRi- fR2{Rl 



i?2 








(A.3) 



where g^ = (g"^ + g'^)/^:. Solutions for the above system are well known [^]. If we 
define t = R2/R1 (obviously t = 1 does not satisfy Eq. ([A.3| )), the solution in the 
ordinary case of radiative SU(2) breaking {m\ + rn^ it 2m| > 0, m\ > m^, rn^m^ — 
m| < 0) is: 



^ {m'f +7712 + 2'm'^){m'f + rrig — 2m 



2m\ 



V 52(t^-l) 
and the global minimum configuration is: 



(A.4) 



i-p arbitrary. 



(A.5) 



f Hi\ ( Ri cos Lp \ 

H2 tRi cos if 

hi Ri sin ip 

\ /i2 / . \ -tRi simp I 

\ / mm \ ^ / 

It is trivial to show that the value at the minimum for this configuration coincides 
with the physical vacuum one of Eq. (|3.6D . 
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B. Tree level constraints for CCB-Eb direction 



Following closely the notation of [11, 12 1, we express all fields in Eq. ( [4.1[1| ) in terms of 
so that 1^1 1 = l-E'il = l3\Hi\, \H2\ = 7|-ffi|- Since the trilinear terms of our example 

have small coupling {\Ye\'^ <C 1), D-terms should be suppressed. This implies that a = P 
and + 7^ = 1. To proceed further, two cases should be examined: 



sign(Ae 



-sign(i?) with mg = fiB, a = y^l — 7^ and < I7I < 1. 



Then, all three trilinear terms can be made negative simultaneously and the tree level 
potential becomes 



(0) 



Ye\'^a^F(a)\Hir - 2YeAa'\Hi\'' + rh'lHi 



c2| 



CCB-Eb 



(B.l) 



where F{a) = 1+2/a^, A = |Ae|+|/i|7 and rh? = m\+{m'\ +m|,Ja^+m|7^— 2|m||7. 
Differentiating with respect to \Hi\ for fixed values of 7, we find besides the trivial 
extremum the following local minimum for Eq. (|B.l| ) 



\H 



SA 



1 ext 



4Yea'^F{a) 



1 + Wi 



8m^F{a) 
9i2 



(B.2) 



Note that the typical vevs are of order |^e|/^e- The corresponding value of the 
potential is 

1 



v. 



ext 



'^,2| TT |2 

l^llezt 



" 21 

m 



YeA\Hi\ext 2" 



(B.3) 



and the constraint to avoid a deeper configuration than the physical vacuum Eq. 
of the theory reads A^ < Fttt? i.e. 

(I^el + 1/^17)' <{l + ^^ {ml + (m|^ + m%y + mW " 2|m^|7) (B.4) 



sign(Ae) = sign(i?) with a = y^l — 7^ and < I7I < 1. 
Similarly the relevant constraints are 



(l^el - < (1 + ^) ("^? + ("^ii + + - 2|"l3l7 



(B.5) 
(B.6) 



C. Field dependent mass matrix elements for UFB-1 

We cite here all the necessary mass matrix elements used in the definition of the one loop 
effective potential. Let ip^ , r = 1,2,3 stands for a Yukawa (Y) or a trilinear soft coupling 
(A). The following notation is used throughout = {ipu, ipc i't) and = {^pd, ips, '4'b)- 
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Gauge Bosons (f2 = C = 1) 
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• Sleptons {L^ = {D i ) n = 2, C = 1) 

'^L = -Y.\^r 1 [Ml]r i; where = ^5 ^2 

r=i \e^J \ai*A'^*aiJ 

We give below the relevant entries 

A[ = Y^i^ihl - Alhi) Al = YJif^H* + AIM,) 

Al = m% + {\H,\' + \h,\') - ^ + \h,\' - \H,\' - M') 

Bl = ml + |y>ip + {\H,f - m') + [M' - M') 

Bl = ml + \Y:H,\^ + ^ - \H,\-) + - M^) 



Squarks (Q^ = ( u J ) O = 2, C = 3) 

3 

{k < ura mi) [M\ 



a,r=l 



( uT \ 



where [M|] 



/ A[ KJ;\ 
Al Kl Kl 
Kl* Kl* Kl K^ 
y K^* Kl* Kl* Kl J 



where 



= + \y:H2\^ + lyj/iip - G- - \H2\^) - G+ {\hr\^ - 



^^2 



Kl 
Kl 



Kl 



9l 



^ - ) nih^ + ( ^ - ii; 



9l 



rr\2 



2ri2 



m 



I + |yj-/fi|2 + |y>2|2 - G+ - |i?2n - G- (|^i|2 - |/i2p) 



Y:;Yj{H,hl + H*2h,) 



mi 



and by definition G+ = {g^ + "igl) /12 G' = {g^ - ?>gl) /12. Besides 



Kl = -Y:;{fiHt + AlH2) 



Kl=YJ{t^H*+A:^H,) 



Higgs (0 = 1, C = 1) 



-^^^[Ml,],^* where 



- 22 - 



with 
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R2 R3 Rq Rj 

214 Uq Ks tig 

, p* p* p* p 
\ -"-5 ^7 ^9 ^10 / 

(^^ ^^*) , = H- Hi) 
R,=ml + ^±^{2\H,\' + \h\'- 



(Qi Q2 Qa Q5 \ 
Q2 Q3 Qe Q7 
Qa Qe Qs Q9 
\ Q5 Qy Q9 Qio J 

Additionally 



R, 



2 2 

R,=ml + ^^{\Hi\^ + 2\hi\'' 


-|^2|'- 


-|/l2p) + ^|i/2p 


Ra = ^H2h\ - ^L_M2LHlh2 


R5 = 


q2 + q2 


p a'^+alr^^u 
Kq — ^ n,x"2 


R7 = 


9^ ~ 92 Tj u* 1 S'2 jLT* I, 


2 2 

p 2 , 5 +52 n U |2 , ol^ |2 

-Kg — m-2 + ^ [\-l^2\ +^1/121 


1 n |2 


1/1 |2\ 1 •S'2 1 77 |2 

- |/iil j + yl^il 








i?io=rn^+^'+^'(2|F2p + |/.2|^ 


-l^il'- 


-|/^il^) + f|/^iP 


and the matrix element for T2 are 








Q2 = 






Qa = 


''-'''H,h2 
4 


Q, = ml ^'^/^H,H2 + ^lh,h2 


Q6 = 


^2 3^ . ,9ljr rr 




Q8 = 




Qg ^ -^2^2 


Qio = 





D. Field dependent mass matrix elements for UFB-3b 

Here we present the mass matrix elements needed in the one loop effective potential 
expression. We use the same notation as in Appendix ^ 

• Gauge Bosons {Q.y/± = 2, ^Ia,z = 1, C = 1) 



1 / A^^ \ 



V2 V3 
V3 Va 



-23- 



i2|' + |i>2|' + |4|') 



'\l3\^-2sm''ew(\El\^ + \i3 



- 1 (1^^ 

2 2 

Vi = [|J/2|^ + |i>2|^ + (l-2sin2eiy)2|4|2 + 4sin^^H^|^||2 

where e = 5(72 / \/ 9"^ + 92 ^-iid sin^^y = 6/52- 

• Quarks (0 = 2, C = 3) 

• Leptons - Higgsinos (0 = 1, C = 1) 

where = (^z/^ r r'^ B j ^nd = (^^r m'' H+ H' W+ W'^ . Besides 



[m(] 



[M 



/ 
















YrEI 


7S'^2 
9 p* 

V2*3 


_92_f.* \ 

92 n* 
^/2^3 













Yrh - 


-V2gEr 











YrE^ 











-A^2 



il^2 






V 


9 r,* 9 f* 
v/2^2 

92 ,-.* 92 0* 


-V29Ef 



V2 2 

92 Tf* 

V2^2 






-Ml 





-M2 






/ 








-YrEl 





-52^1 ^ 


















-52 ^^2 


















-Y^i'2 





















-M 


- 


92H*^ 






-Yrh 


-Yf,U2 


-M 















-92V*2 











-M2 






V -92tl 





-52^2* 





-M2 








Squarks {CT = { u d ), = 2, C = 3) 



= -E 

r,a=l 



r,o=l 



f 1 and [M|] 



where [Mf] = 

if[ = + |y:iJ2|' 




4r 



. The relevant quantities are 

^-^{2\Elf + \H,f-\i>,f-\isf) 

9^ - ^9l 
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• Sleptons - Higgses {L^ = ^ )) 



where y-^ 



M2 = 



= i>3 H-* ii) and </,^ = {e e*) , e = HI El . 

4, + y (l^~3p - \Elf) + ^ (|^2p - 141^ - |i^2p) (J^ = 2, c 

On the other hand [M^ ] = [ | (j^ = 2, C = 1) 



' = 1) with 



A. = 



m 



a_ 

2 

1 + ^ (141^ - l^af ) + ^ (l^2|^ - |4P - |i^2p) . 

/ 5i B2 B3 -ml B4 \ 
B2 B5 Bq B-j Bg 

rj* D D D 

-03 -De -°9 -DlO -Dll 

-mi S| B*Q Bu 
V S| B* B*^ 

Si = mi + ^ {\Elf - + {ihl' + |J/2p - 



Similarly [M|J = 



(0 = 2, C = 1) with 



B^ = ml 



Be = Y^YrU2Ef Bj = -Y^A^v^ 



Bs = Yi^ifiH^ + YJlE^) 



Bo 



c,2 , 9 +52 ^,^^|2^|4|2_|^^,2 



(I 



) 



BiQ = —Y-tAt-E' 



Bn = 4^*21*3 



Bu = ml + |y,P2p + \YrEI\' + Y (l^~3p - + {ihl' - |4P - \H2^' 

2 , I rc|2^ 



Bi3 = ml + 



z 2 

2 I g + g2 A--, |2 I |« |2 I I Tj |2 
Z2 + ^ [l^^l + KSI + l^2| 



EL 
2 



-f^2 + -E'fl 
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and [M 



S3 J 



^1 ^2 I (^Q^c=l) with 



f Ri \ 

i?2 ^3 ^4 ^5 
i?3 i?6 ^7 ^8 
i?4 -R7 Rg RiQ 

\0 RlRl Rio RnJ 



Do 



/ m| Qi Q2 \ 
ml Q3 Qi Q5 Qe 
Qi Qi Qi Qs Q9 
Q5 Qs Qw Qu 

\ Q2 Qe Q9 Qu Q12 1 



Ei = m? + |y.p(|^3f + l4P)+|-(|^3 

-.2 



i?2 

Rq 

R7 



R 



10 



\'-\m\') + ^—^(\^2\' 



|2\ , 5^+52 



+ 



4 

|z^2| 



Ra 



|J^2 
^2 



Ro = mi + 



|2 _ |^.|2^ + 9^ j^2|^>2p - |4P - |^2P) 



L2 



^ ^2^3 



1^1' - l^2|' 



2^3 



^2P 



2^31' - 



while the corresponding elements for D2 are 

Q2 = Y^ArEl* 



■,2 



Q4 



a 



-^2^3 



= H2l'2 



Q^='^^{H2? 

Q, = -l^hH2 



Q7= (gElf 

Q9 



^,YrH2+{\Yr\^-^]i3El 



Q\\ = — z — 1^213 



Qs-- 

Qio 

Q12 



g"^ +92r~ a2 



4 

^2 I „2 



{^2y 



E. Field dependent mass matrix elements for CCB-Eb 



Using the notation of Appendix the mass matrix elements needed in the one loop 
effective potential expression are 
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Gauge Bosons {0,\y± = 2, flA,z = 1, C = 1) 



V2 



9l 



141 



C|2 



m2 = |y(.ii-i| 



where e = 552/ \J 9^ 9\ and sin^^ = 6/52- 

• Quarks (J^ = 2, C = 3) 

• Heavy Leptons (fi = 2, C = 1) 

• Light Leptons - Higgsinos (0 = 1, C = 1) 

£2 = -^VfKlV'l " ^^2'^[M/]V2 

where ^\ = (e e" B 1^(3)^ ^ ^^.^ ^- ^+ ^+ ^^g^ 



[Ml. 



( ° 









S /* 32 /* 
V2*l V2^1 












-V25^r 




ye^i 







9 TI* 92 TJ* 

73^1 73^1 













9 TJ* 92 TJ* 


9 p* 


-\/2gEl* 


_2_ rr* 




-Ml 


92 f* 





92 TJ* 


71^2 


-M2 



[Mi] = 



( -ye^f -g^ll \ 

-Y^El -II -g2H* 
-At -52i^| 

-ff2i?i* -M2 

V -52^1 -52i^| -M2 y 



Squarks (Q^ = ( u d), n = 2, C = 3) 



= - E (f^™ [^1] ( ^? ) - E [Mi] 



r,a=l 



r,a=l 



mi 
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where [Mf ] = J. and [M|] = J I The relevant quantities are 



r" \ A^* A^ 



Tl = + \Y:;H,\^ + ^ - U) + (^R,f + U - 

AS = + iFjiJip + ^ (l^^p - + [\H,\^ + - \H,\^) 

• Sleptons - Higgses {L^ = ^ )) 

= -6l^2p - 6I^3P - (^i ii) [Ml] I* ^ 

- ii) [Ml] l^f^- x^[Ml]x* - \f[Ml]ct,* 

where = {H^ H+* i>*) and cf^ = ((^^ ^t) ^ ^ j^^o ^0 
For p = 2, 3 using a compact notation we have 

= + ^ - + (iifip - - {n = 2,c = i) 

On the other hand [Ml] = l^l j = 2, C = 1) with 



Z2 = y;,(/xif2* + A^Hi + YJiEi*) 



Z, = m|. + \Y,H,\' + ^ (|^2|' + 2|^^|2 - 141^ - \H,\') 



Z, = ml + \Y,H,\' + ^- - I^Jp) + (lifip - |i72p - \£i\') 



^^d[Ml]=l^'^^'] (Q = 2, C = 1) with 



K2 = yr(Ai^2* + AHi + YeilEl*) 



Ks = m^- + + ^ f _ |^C|2\ + , |^^|2 _ |^^|2 _ 
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Bi B2 -B3 \ 
Similarly [M^] = | 5* ^4 

B*3 Bl Be ) 



(Jl = 2, C = 1) with 



Si 
B4 



mo + — 
2 2 



5fi 



< + (l^ep-|-)|^lf + 



5^ A I i^c\2 , + fl'i / I |2 , I rr |2 I tr |2 

1 + -ni — Ln2 



and [M2] 



^2 SI 



{n = C = 1) with 



/ 61 62 63 64 \ 
65 96 67 

Gg Gg 

VG| G^ G^ Gio/ 



So = 



/ Ai A2 A3 A4 \ 
A2 A5 Ae A7 

A3 As Ag Ag 
\ A4 A7 Ag Alo ) 



Gi = m? + lYeP {\El\' + \h\') + \ {^h? - \E\f) + (: 



62 = -l±AHtH2 



e^={\Ye\''-^-]HlEl 

9^ + 92 



Q.= \\YA' + '—r^ 



Ge = M^e^t + \H*El 



G7 = /xFe^r - 



c* 5^ - 5i 



2*1 



©9= |yeP-^U~i^? 



Gio = m|^ + {\HA' + |^^|^) + ^ [2\h\' - \El\') + - |if2|' - 2|4|') 

while the corresponding elements for ^2 are 



A, = ^^[H,f 



A2 = ml-^±^H^H2 



A3= {\Y,\^-^]HiEt + Y,AXi 



A4=Ye(YeHji+AeEr) 



9^-9, ~ 
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At 



2^1 



Ag = Ye inH2 + AeHl) + ( \Ye\^ - ^ ) hE^ 



A6 = 
A8 = 

Alo 



9 



2 



HoE' 



2^1 



I? 



F. Expressions for uj{x) 

We present below expressions of u){x) used in the definition of the renormalization 
scale Q* given in Eq. (|5.4D for the various cases we have considered 



• UFB-1 



ijj{x) 



All (^12 + x)(Ai3 + x)(Ai4 + X) 
^21(^22 + X){A2Z + X){A2A + x) 

Azi{A:i2 + x)(A33 + x){A^i + x) 

^41(^42 + 2;)(A43 + ^442; + x"^) 
^51 (^52 + X){A^^ + A^iX + 
^61(^62 + a;)(A63 + ^643; + X^) 
^71 (^72 + x){Aj^ + ^743; + 
^81(^82 + 2;)(A83 + ^842; + X^) 
^l(^2 + X)(^3 + ^94X + 



where 



• UFB-3b 



uj{x) 



X < 2.8 
2.8 < X < 3.0 
3.0 <x< 4.0 
4.0 < X < 5.0 
5.0 < X < 6.0 
6.0 < X < 7.0 
7.0 < X < 8.0 
8.0 < X < 9.0 
9.0 < X 



/ -8.22877 


-3.07458 


-2.70973 


-2.3097 \ 


1.40049 


-3.96386 


-3.51934 


-2.71471 


0.86922 


-4.78838 


-3.43647 


-2.70523 


-0.390767 


-3.5568 


30.1795 


-10.8231 


0.0938483 


-8.08594 


24.5516 


-9.49597 


-0.184626 


-4.67762 


48.0765 


-13.5349 


0.344656 


-9.61596 


42.2498 


-12.8772 


-0.193997 


-6.81241 


102.729 


-19.8646 


V 0.0313323 


-14.349 


63.3896 


-14.651 / 



r Uu{ui2 

f/2l(f/22 
f/3l(f/32 

Ujl{Uj2 
Usi{U82 
Ugi{Ug2 



+ X)([/13 
+ x)([/23 
+ x)([/33 
+ X)(?743 
+ x)(C/53 
+ X)(?763 

+ x){Ur3 

+ X){U83 
+ x){Uq3 



+ x)([/i4 + x) 
+ x)(f/24 + x) 
+ x)([/34 + x) 
+ [/44X + X^) 
+ f754X + X^) 
+ f764X + X^) 
+ f774X + x2) 
+ f784X + X^) 
+ f794X + X^) 



X < 2.8 

2.8 < X < 3.0 

3.0 < X < 4.0 

4.0 < X < 5.0 

5.0 < X < 6.0 

6.0 < X < 7.0 

7.0 < X < 8.0 

8.0 < X < 9.0 

9.0 < X 



(F.l) 



(F.2) 



(F.3) 
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where 



U 



• CCB-Eb 



uj{x) 



where 



C 





_Q 1 1 QCiq 
O . -L X ty 




— 9 9fi'^n9 \ 


— n Qnni 






— 1 9049^ 


U.VJOtJ 1 '-tU 




— 3 39069 


— 9 81 41 8 


-0.292473 


-3.53026 


28.8914 


-10.313 


0.184148 


-8.25836 


19.5607 


-8.57881 


-0.144118 


-4.62537 


61.7438 


-14.8604 


0.192324 


-10.2975 


37.4871 


-11.8372 


-0.125179 


-6.15239 


114.679 


-20.7135 


V 0.00839312 


-18.1679 


59.3748 


-10.8321 / 



Cll(Cl2 
C2l{C22 

C4i(C42 
C5l(C'52 
C6i{Cq2 
C7l{Ci2 
Csi{Cs2 
, C'9l(C'92 



+ x)(Ci3 
+ X)(C23 
+ X)(C33 
+ X)(C43 
+ X)(C53 
+ x){CQ:i 
+ X)(C73 

+ a::)(C83 
+ a;)(C93 



+ x)(Ci4 + x) 
+ C2iX + x^) 
+ X)(C34 +X) 
+ C44X + x^) 
+ C54X + 

+ C64a; + x^) 

+ C74X + 

+ (784^ + X^) 

+ C94X + x^) 



X < 2.8 

2.8 < X < 3.0 

3.0 < X < 4.0 

4.0 < X < 5.0 

5.0 < X < 6.0 

6.0 < X < 7.0 

7.0 < X < 8.0 

8.0 < X < 9.0 

9.0 < X 





-8.31888 


-3.07077 


-2.70937 


-2.31018 \ 




1.9263 


-2.71468 


12.4319 


-7.02268 




0.581498 


-5.27073 


-3.4845 


-2.68738 




-0.0667458 


-3.55197 


53.9364 


-13.3041 




-0.0145147 


-3.39174 


127.123 


-20.1434 




-0.0751954 


-3.95696 


65.0836 


-15.1115 




0.0152963 


-19.2871 


39.5677 


-11.2084 




0.11401 


-11.7758 


46.3403 


-13.0956 


v 


-0.0713367 


-6.7888 


150.536 


-23.613 / 
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